Draft version December 26, 2009 

Preprint typeset using I4TgX style emulateapj v. 08/22/09 



IDENTIFICATION AND PROPERTIES OF THE PHOTOSPHERIC EMISSION IN GRB090902B 

F. Ryde 1,2 , M. Axelsson 3,2 , B. B. Zhang 4 , S. McGlynn 12 , A. Pe'er 5 , C. Lundman 2,1 , S. Larsson 6 ' 2 , M. Battelino 1 



B. Zhang 4 , E. Bissaldi 7 , J. Bregeon 8 , 

F. LONGO 1314 , S. MCBREEN 7 ' 15 , 



, M. S. Briggs 9 , J. Chiang 10 , F. de Palma' 1 ' 2 , 
N. Omodei 12 , V. Petrosian 16 , R. Preece 9 

Draft version December 26, 2009 



S. Guiriec 9 , J. Larsson 32 , 



A. J. VAN DER HORST 



17 



o 
o 

o 

Q 

(N 

43 

9 L1, 
6 



> 
in 

o 



o 



ABSTRACT 

The Fermi Gamma-ray Space Telescope observed the bright and long GRB090902B, lying at a redshift of 
z = 1.822. Together the Large Area Telescope (LAT) and the Gamma-ray Burst Monitor (GBM) cover the 
spectral range from 8 keV to >300 GeV. Here we show that the prompt burst spectrum is consistent with 
emission from the jet photosphere combined with nonthermal emission described by a single powerlaw with 
photon index -1.9. The photosphere gives rise to a strong quasi-blackbody spectrum which is somewhat broader 
than a single Planck function and has a characteristic temperature of ~ 290 keV. We model the photospheric 
emission with a multicolor blackbody and its shape indicates that the photospheric radius increases at higher 
latitudes. We derive the averaged photospheric radius i? p h = (1.1 ± 0.3) x 10 12 F 1 ' 4 cm and the bulk Lorentz 
factor of the flow, which is found to vary by a factor of two and has a maximal value of T = 750 Y 1 / 4 . Here 
Y is the ratio between the total fireball energy and the energy emitted in the gamma-rays. We find that during 
the first quarter of the prompt phase the photospheric emission dominates, which explains the delayed onset 
of the observed flux in the LAT compared to the GBM. We interpret the broad band emission as synchrotron 
emission at R ~ 4 x 10 15 cm. Our analysis emphasize the importance of having high temporal resolution when 
performing spectral analysis on GRBs, since there is strong spectral evolution. 

Subject headings: gamma rays: bursts - gamma rays: observations - gamma rays: theory - radiation mecha- 
nism: thermal 



1. INTRODUCTION 

The mechanism giving rise to prompt emission in gamma- 
ray bursts (GRBs) has long been a puzzle. The emission 
typically peaks in the 100-1000 keV range and is mod- 
eled with the empirical Band functio n, which consis ts of 
two exponentially joined power-laws (Ban d et al.lll993l) . In 
general the narrow spectral ranges available during obser- 
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vations of GRBs have made it difficult to unambiguously 
determine the emission mechanism. This has been reme- 
died by the launch of the Fermi Gamma-ray Space Tele- 
scope which regularly observes GRBs with its two instru- 
ments, the Large Ar ea Telescope (LAT, n ominal energy range 
20 Me V->300 GeV; lAtwood et aU2009[) and the Gamma-ra y 
Burst Monitor (GBM, 8 keV-40MeV: iMeegan etaH 120091) . 
thus covering an unprecedented spectral range. 

The bright, long GRB090902B was detected by Fermi 
and was one of th e brightest observed to date by the LAT 
dAbdo etalJ 120091) . Over 200 photons were detected at en- 
ergies above lOOM eV, including 39 photons with energy 
above 1 GeV dde Palma et alj |2009). The LAT also detected 
a photon with energy ZZA^^ GeV, the highest seen in any 
GRB so far. The bu rst lies at a redshift of z = 1.822 
(Cucchiara et al.l 120091) . which yields an isotropic energy of 
£-iso = (3.63 ± 0.05) x 10 54 ergs, and the afterglow emission 
was seen at X-ray, optical, NIR and radio wavelengths. 

The prompt spectrum of GRB090902B over the energy 
range 8 keV-33 GeV shows c lear deviation fro m the gener- 
ally expected Band function ( Abd o et all 120091) . The time- 
integrated spectrum is best fitted by the addition of a sepa- 
rate power-law component (photon index of ~ —1.9) to the 
Band function. The power-law is detected at both lower and 
higher energies around the Band component. During the first 
half of the prompt emission phase, the Band power-law index 
at energies below the spectral peak significantly violates the 
optically thin synchrotron limit; photo n index of a = —2/3 
dPreece et al. Ill 9981; iPreece et alJEOOH) . 

Rvde (2005) suggested that GRB spectra are a superpo- 
sition of two spectral components: photospheric blackbody 
emission and an accompanying nonthermal component. Over 
the limited energy range of 20-2000 keV (using data from 
BATSE on CGRO) studied, an adequate model consists of a 
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single Planck function and a power-law. However, the non- 
thermal component, which is modelled by the power-law, is 
not well characterized, since the burst spectral shape is i n gen- 
eral dominated by the therm al component (see further Rvde 
(120081) : iRyde & Pe'eri (l2009h V The Fermi data now aUow an 
improved characterization of the two components due to the 
increased energy range. In this paper, we investigate in more 
detail whether the observed two components in the broadband 
spectrum of GRB090902B could be attributed to such a pho- 
tospheric model. 

2. TIME-RESOLVED SPECTRAL ANALYSIS 



lAbdo et af] ((2009) analysed time-resolved spectra of 
GRB090902B with time bins of typical duration of 6 s. In 
particular, they studied the behavior of the power-law compo- 
nent in the spectrum. In this work, we investigate the peaked 
component of the spectrum, i.e. the "Band" component, at a 
higher time resolution in order to study the spectral evolution 
in greater detail. A consequence of this i s that any model w ill 
be less constrained at LAT energies (cp. lAbdo et al.1 ((2009)). 
Of particular interest is the first half of the prompt phase 
(t = — 13 s), during which the Band component is very hard 
and narrow, which challenges optically-thin emission mech- 
anisms as its origin. The Band component has a low-energy 
power-law Ne oc E a — E° and a high-energy power-law 
N E oc E 13 = E~ 4 . We theref ore focus our study on epochs 
a, b, and c in lAbdo et al.l (120091) . and only briefly comment on 
the later behavior. 

The time binning was chosen by requiring a signal-to-noise 
ratio of 40 in the most strongly illuminated GBM detector, 
Nal 1 (see Table 1). We also include data from Nal detector 
0, BGO and 1, and from the LAT in our fits. The LAT 
"transient" class data c ontain front and bac k events which are 
considered separately (lAtwood et al.ll2009l) . The Nal data are 
fit from 8 keV to 1 MeV and the BGO from 250 keV to 40 
MeV using the Time Tagged Event (TTE) data type. The LAT 
data are fit from 100 MeV. An effective area correction factor 
of 0.9 is applied to the BGO with respect to the Nal detectors 
and LAT. The fits were performed with the spectral analysis 
software package RMFIT (version 3.0). 

The first result we find is that the a-value gets larger, i.e. the 
sub-peak spectrum gets harder when narrower time intervals 
are used fo r the Band+power law model. For the time inter- 
vals used in lAbdo et al.l (12009) a ~ 0; with the narrower time 
bins adopted here, we obtain an average value of a ~ +0.11 
(see Table 1). Several spectra have a ~ +0. 3, which are 
among the hardest GRB spectra ever measured (Ryda l2004t 
iKaneko et alJl2006|) . We therefore attempt to interpret this 
hard and narrow spectral component as stemming from the 
photosphere. First we model it with a single Planck func- 
tion, which is used in addition to the powerlaw (BB+pl). 
The average value of C — stat/doi = 634/599 (reduced 
C — stat = 1.058). We note however that the distribution 
of the residuals between the model and data shows trends that 
are not expected from stochastic variations, and thus indicate 
that the peak is slightly broader than a single Planck function. 

Indeed, from a theoretical point of view the photospheric 
emission is not expected to be a pure Planck function. A 
broadening of the thermal component is expected due to con- 
tributions from different regions in space. Goodman (1986) 
discussed the emerging spectrum from an optically thick, 
fully ionized relativistic flow and showed that the emerging 
spectrum has a peak that is broader and a slope below the 
peak that is slightly shallower compared to the Planck func- 



tion. Several effects needs to be taken into account. The ob- 
served blackbody temperature depends on the latitude angle 
due to the angle dependence of the Doppler shift. Likewise, 
the optical depth is angle dependent, which results in the pho- 
tospheric radius increasing with angle (Pe'er 2008). Angle 
dependent density profiles of the outflow will have a similar 
effect. Therefore, the photospheric component is better rep- 
resented by a multicolor blackbody instead of a single Planck 
function. The broad energy range now available through the 
Fermi observations, allows us to model the spectrum of the 
photospheric emission in greater detail than previously since 
the accompanying nonthermal component can now better be 
constrained. We thus fit the photospheric component with 
a phenomenological multicolor blackbody (mBB), which is 
given by 



F mBB (E,T u 



E A 



dA(T) 
dT exp[E/kT] - 1 



dT 



(1) 

The spectrum thus consists of a superposition of Planck func- 
tions in the temperature range T = T m [ n to T max ; T max is 



a free parameter (T n 



« T„ 



and therefore cannot be 



determined). The spectrally integrated flux for each Planck 
function is given by F(T) = A(T) T 4 tt 4 /15, where A(T) is 
the normalization. Phenomenologically, we introduce the in- 
dex q, relating the flux and the temperature of the individual 
Planck functions 



F(T) = F„ 



T 



(2) 



where F max = F(T — T max ). This full thermal + nonther- 
mal model has 5 free parameters; power-law index, powerlaw 
normalization, T max , ^4(T max ), and finally q. 

Figure [TJ shows the time-resolved E F-& spectrum from one 
of the time intervals (1 1.008-1 1.392 s after the GBM trig- 
ger). The mBB function captures the sharp spectral peak 
emerging above the nonthermal component, being a factor 
of ~ 30 above it. Such a strong therm al component is 
rarely seen in GRB spectra (Ryda [20041) . For this time- 
bin feT max = 270.lt\z l keV, and K = (F hh /aT 4 ) 1/2 = 
(3.50 ± 0.05) x 10" 19 dPe'er et al.ll2"007h . The power-law 
index of the nonthermal component is — 1.951q'oi3 and the 
g-parameter 2.00^q'^2- The fit has a C — stat/doi = 
516/598 = 0.86. The quality of the fits improve com- 
pared to the BB+pl model. Over all time bins the averaged 
C — stat/doi — 0.944, which is an improvement of C-stat 
with 69 for 598 dof. Moreover, the residuals do not have any 
marked trends. 

In Fig. |2] the evolution of the characteristic temperature is 
shown. Due to the brightness of the burst and the relative 
strength of the thermal component, the temperature is deter- 
mined with unprecedented accuracy. Since the burst consists 
of heavily overlapping pu lses the characterist ic cooling be- 
havior over a single pulse dRvde & Pe'erll200"9l) is not clearly 
detected. The middle panel in Fig. [2] shows the dimension- 
less parameter 7Z. For the studied interval in GRB090902B 
it is fairly constant, varying only b y a factor of 2. Typ ically, 
1Z varies by an order of magnitude (Ryde & Pe'er 2009). Fi- 
nally, the g-parameter does not vary much during the analyzed 
interval and has an averaged value q = 1.9 with a standard de- 
viation of 0.4. 

An alternative interpretation for the broadening of the ther- 
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mal component is the effects of scattering and Comptoniza- 
tion, as worked out in IPe'er et al.1 (120051120061) . Moreover, 
heating of thermal electrons by Coulomb collisions with pro- 
tons will give rise to a broadened photospheric component 
dBeloborod ov 2009). We therefore also fit the photospheric 
component with a spectrum taking these effects into account. 
This is approximated by a Band function with a = 0.4 and 
/3 = —2.5. The averaged C — stat /dof = 0.937, is similar 
to the mBB model and the residuals do not have any marked 
trends. 

The nonthermal component is adequately fitted by a single 
power-law and a curvature in the spectrum is not statistically 
required. The photon index is constant, ~ —1.95, with a stan- 
dard deviation of 0.05. This value is largely independent of 
the width of the ti me bins as well as of model used for the 
thermal component (Ab do et alJl2009h . 

The fraction of energy flux in the thermal emission relative 
to the total flux in the observed energy band is on average 
70%. During the first 6 s the average is 86% (the spectrum 
is dominated by the thermal comp onent), after which it set- 
tles to an averag ed value of 63%. Abd o et al.l d2009l) noted 
that there is a delayed onset of the LAT light curve compared 
to the GBM light curve. This can be naturally explained by 
the dominance of the photospheric component during the first 
6 s of the burst. To investigate this further we reduced the 
signal to noise ratio to 10, increasing the number of time in- 
tervals. The q-parameter was frozen to 1.9 (the average value) 
and the data were fit with the photospheric component, with 
only two free parameters, the temperature and the normaliza- 
tion. The resulting, averaged reduced C — stat has an accept- 
able value of 1.02 for t = — 6 s. Thereafter the reduced 
C — siai-value increases dramatically to ~ 4 , indicating 
that an extra component is needed, i.e. the power-law com- 
ponent. This is thus a natural explanation for the observed 
lack of photons in the LAT and in the Nal below 14 keV at 
the beginning of the pulse. We note that a similar behav- 
ior was observed by BATSE i n GRB 970111 (Fig.13 in lRvdd 
(2004T) lGhirlanda et all (120041) ). 

Finally, we calculated the cross correlation function be- 
tween the thermal and nonthermal light curves, which are 
found by integrating the mBB and power-law functions over 
the observed energy band for every time bin. We find a lag 
of 0.47 ± 0.27 s (thermal emission leading) and a correlation 
maximum of 0.91 ± 0.09, indicating a strong correlation. 

Beyond t he t ime interval studied here, we showed in 
lAbdo et aTl (120091) that the Band component becomes broader 
and softer. We find that the photospheric model (mBB) can 
adequately fit these spectra as well, however the value of the 
q parameter is markedly lower, lying between 1 and 1.5: the 
properties of the photosphere changes. Further theoretical 
study of the interpretation and reason for this is underway. 

We thus conclude that by resolving the light curve on a sub- 
second timescale one can identify a photospheric component 
in the spectrum, combined with a single power-law compo- 
nent. The photospheric component can be modelled by a hard 
and narrow Band function. A physical interpretation is given 
by a multicolor blackbody, which has on e parameter less t o 
be fit. Using wider time bins, as done in ( Abdo et al.ll2009h . 
will increase the SNR but the variations in temperature and 
flux will broaden the spectrum, weakening the thermal signa- 
ture. Our results thus emphasize the importance of having the 
highest possible temporal resolution when performing spec- 
tral analysis on GRBs. 



2.1. Properties of the photosphere 

The identification of the emission from the photosphere al- 
lows us to determine physical properties of the relativistic out- 
flow, such as the bulk Lorentz factor T, the photospheric ra- 
di us i? D h and the ini tial size of the flow Rq, as we showed 
in IPe'er et ail d2007l) . The photospheric radius is given by 
Rph — Lo-T/8irT 3 m p c 3 , where L is the total fireball lumi- 
nosity, and <tt and m p are the Thomson cross section and 
proton mass respectively. The values of i? p h and T can be 
determined by combining this expression with the measure- 
ment of the dimensionless parameter 1Z, which is related to 
the effective transverse size of the photosphere 

K={ <I±i>!3* (3) 

d L r 

where £ is a geometrical factor of order unity, and is the 
luminosity distance. Furthermore, we can estimate the radius 
above which the relativistic acceleration begins, Ro. This 
is done by combining classical, nondissipative fireball dy- 
namics with the fact that the total fireball luminosity L = 
And\Y F^s, where Y is the ratio between the total fireball 
energy and the energy emitted in the gamma-rays, where 

-^obs — -^bb ~t~ -^non-th ■ 

The discussion in IPe'er et aTl (120071) assumed a single 
Planck spectrum. However the discussion can be generalized 
to slightly distorted photospheric spectra, as in the ones dis- 
cussed above. For instance, the q-value we find here (q ~ 2) 
indicates that at a certain observer time the photospheric ra- 
dius increases with latitude angle, so for every time bin we 
therefore estimate an average photospheric radius. The values 
found below are similar to the ones found by approximating 
the spectrum with a single Planck function and applying the 
theory directly. 

All the parameters evolve with time. In particular we 
find that the bulk Lorentz factor T starts off at a value of 
close to r ~ 550 for the first 6 s. Therea fter it rises 
sharp ly to a maximal value of T ~ 750 (cf. Abd o et al.l 
2009). The time -averaged value (with standard deviation) 

for the bulk Lorentz factor is Y = 580 ± 130(£F) 1/4 , 
while the time-averaged photospheric radius is found to be 
i? ph = (1.1 ± 0.3) x lO^cm^ 3 / 4 ^^ see the right panel 
in Figure [2] The photospheric radius is remarkably stable 
compared to the variations in temperature and flux. Further- 
more, since T is most strongly dependent on the tempera- 
ture, the evolution of these quantities track each other. We 
also estimate the radius at which the jet is launched to Rq — 
(1.0 ± 0.5) x 10 9 £~ 4 y~ 3 / 2 cm, and consequently the satura- 
tion radius R s = TR Q = (5.2±1.7) x 10 11 £- 15 / 4 y- 5 / 4 cm. 
Here we assume nondissipative fireball dynamics. However, 
we note that the dynamics d epends on the magnetization 
dGiannios & Spitkovsky| [2009) and on the presence of sub- 
photospheric heating, such as tangential col limation shocks 
dRees & Mesza rosll2005l lLazzati et al. 1120091) . and Coulomb 
heating (Belob orodovl2009l) . 

3. DISCUSSION 

The presence of a strong photospheric component in GRB 
spectra has been discussed by several authors (e.g. iGoo dman 
[l98llMeszaros et al. ll2002l;lRees & Meszarosll2005l) 18 . Here 

18 se e also e.g. IPaczvnskil 119861): Meszaros & Rees 

(2003); ILvutikov & Usovl (2000); Daigne & Mochkovitch 120021) : 
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we have identified a strong photospheric component in 
GRB090902B, which is modelled by a multicolor blackbody. 
It dominates over the nonthermal component, being nearly 
100% at early times. 

Varying energy injection at the central engine, or variations 
caused by interaction with the progenitor material can cause 
the observed variability in the light curve. The flow is ad- 
vected through the photosphere, where the thermal emission 
escapes. A fraction of the kinetic energy stored in the flow 
is later dissipated and emitted as nonthermal emission, e.g. 
synchrotron radiation. It is therefore expected that both the 
thermal and nonthermal emission reflect the properties and 
original energy injection at the central engine and thus are 
correlated with each other. This is consistent with the high 
correlation observed, with a time lag of approximately 0.5 
s. Such a lag corresponds to a shock radius i? s h — R p h = 
2cr 2 i lag (l + zY x - 4 x 10 15 (r/580) 2 cm ~ R sh at which 
the optical depth r ~ (-R s h/-R P h)~ 2 ~ 10~ 7 ensuring an op- 
ti cally thin-emission site. 

Abd o et aT] (|2009) estimated a lower limit on the bulk flow 
Lorentz factor r m j n ~ 1000 by combining the energy of the 
11.16 GeV photon, observed at 11 s after the GBM trigger, 
with the variability timescale in the LAT data which was de- 
termined to be t v ~ 0.1 s: The gamma-ray opacity for pair 
production should be less than unity in order to allow the 
photons to escape without attenuation from the flow. This 
Lorentz factor is slightly higher than the value we determined 
above for the flow advected through the photosphere. How- 
ever, the estimations of T based on opacity arguments are 
sensitive to uncertain ties in determining the variability time 
(Zhan g& Pe'erll2~009t) . Indeed, assuming that the variabil- 
ity timescale in the nonthermal component is set by the an- 
gular timescale at the shock radius, we have t v = = 
(1 + z) i? s h/(2cT 2 ) = 0.5 s. This value is consistent with the 
fact that there is not much variability strength on timescales 
shorter than 1 s, as measured by the power density spec- 
trum of the LAT data. A variability timescale of 0.5 s yields 
r m in ~ 750, which is close the value we estimate. Finally, 
we note that the observed variability timescale in the GBM 
data should reflect the angular timescale of the photosphere. 
The determined value tf BM ~ 50 ms (lAbdo et all 120091) is 
larger than (1 + z) i? ph /(2 cT 2 ) ~2x 10~ 4 s, which is to be 
expected as shown by Pe'er (2008). 

As shown above, the photospheric emission allowed us 
to estimate Rq assuming nondissipative acceleration: Ro ~ 
10 9 £ _4 Y _3 ' /2 cm. This is larger than the physical size around 
a solar mass black hole, typically assumed to be at the 
centre of the GRB engine , with a Schwarz schild radius of 
R sch = M /c 2 ~ 10 7 cm (lPaczynskilll986T) . This might be 
an indication of that the fireball dynamics includes signifi- 
cant dissipation during the acceleration phase and/or subpho- 
tospheric heating, which would lower the estimated value of 
i?o- Alternatively, the larger size we find might correspond 
to the radius of t he stel lar core, at which the jet is launched. 
Thompson et al.l (120071) argue that internal dissipation within 
the star prevents the Lorentz factor from significantly rising 
until the jet escapes the co re of the progenito r star, such as a 
Wolf-Rayet star. However. lLazzati et al. I (120091) ) showed that 
the dissipative shocks in the jet continue far beyond the stellar 
radius, thereby influencing the jet dynamics. 



Some bursts have no strong evidence of a narrow thermal 
component. For cases like GRB080916C which has a stan- 
dard Band spectrum covering 6-7 orders of magnitude, the 
lack of a thermal component can be taken as an argum ent for 
a Poynting flux dominated flow (Zhan g & Pe'erl l2009). For 
other cases that show a standard Band spectrum in a narrower 
energy range, like during the second half of the prompt phase 
in 090902b, it is possible that a modified photosphere spec- 
trum with subpho t ospheric heating ca n account for the dat a 
dPe'er et al.l 120071: iBeloborodovl 120091: lLazzati et al.l 120091) . 
Broadband observations of more GRB prompt emission spec- 
tra are desirable to more definitely diagnose the composition 
of GRB jets. 
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FIG. 1. — Time-resolved vF v spectrum for the interval t = 11.008 — 11.392 s over the GBM + LAT energy ranges fitted with a multicolor blackbody + 
power-law model. 
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FIG. 2. — Evolution of (a) kT, (b) TZ, (c) V (black) and i? p h (red). The yellow line indicates the time of the detection of the 1 1.21 GeV photon. 
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Spectral fits. 
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